Two-fluid and population balance models for subcooled boiling flow

Population balance equations combined with a three-dimensional two-fluid model are employed to predict subcooled boiling flow at low pressure in a vertical annular channel. The MUSIG (MUltiple-SIze-Group) model implemented in the computer code CFX4.4 is further developed to accommodate the wall nucleation at the heated wall and condensation in the subcooled boiling regime. Comparison of model predictions against local measurements is made for the void fraction, bubble Sauter mean diameter and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Additional comparison using empirical relationships for the active nucelation site density and local bubble diameter is also investigated. Good agreement is achieved with the local radial void fraction, bubble Sauter diameter and liquid velocity profiles against measurements. However, significant weakness of the model is evidenced in the prediction of the vapour velocity. Work is in progress to circumvent the deficiency through the consideration of additional momentum equations or developing an algebraic slip model to account for bubble separation.     

Study on the possibilities of controlling the laser output beam properties by an intracavity deformable mirror

The possibilities of controlling the laser beam properties by a deformable mirror introduced into the laser optical cavity were studied theoretically and experimentally. The experiments were performed under conditions of an industrial high power transverse flow cw CO₂ laser operating with a stable resonator of a folded configuration. A deformable bimorph mirror of a surface profile controlled by the voltage applied to the mirror electrodes is implemented to the laser system as a back cavity mirror or as a one of the inner folding mirrors. The near-and far-field characteristics of the laser beam versus the resonator configuration controlled by the changes of the focal length of the deformable mirror are discussed in the paper. The analysis reveals that the resonator with an inner deformable mirror is much more sensitive to the mirror curvature variations than the resonator in which the deformable mirror is used as a back cavity mirror. The presented results show that dynamic and controllable changes in the resonator properties result in the controlled modification and optimisation of the laser output power and spatial parameters of the laser radiation.     

The flow of a variable viscosity fluid between parallel plates with shear heating

A model for the flow of a fluid through a channel with parallel plates is investigated. The channel is narrow, so that the lubrication approximation may be applied. The channel walls are maintained at a constant temperature. Shear heating effects are included and the fluid viscosity decreases exponentially with temperature. When the flow is driven solely by shear stress or imposed velocity at the top, analytical progress is possible. When pressure gradient also drives the flow the problem is solved numerically.     

Integration and propagation of a multi-criteria decision making model in constraint programming

In this paper we propose a general integration scheme for a Multi-Criteria Decision Making model of the Multi-Attribute Utility Theory in Constraint Programming. We introduce the Choquet integral as a general aggregation function for multi-criteria optimization problems and define the Choquet global constraint that propagates this function during the Branch-and-Bound search. Finally the benefits of the propagation of the Choquet constraint are evaluated on the examination timetabling problem.     

Dynamics of molybdenum nano structure formation on the TiO2(1 1 0) surface: A kinetic Monte Carlo approach

The rutile TiO₂(1 1 0) surface is a highly anisotropic surface exhibiting “channels” delimited by oxygen rows. In previous experimental and theoretical DFT works we could identify the molybdenum adsorption sites. They are located inside the channels. Moreover, experimental studies have shown that during subsequent annealing after deposition, special molybdenum nano structures can be formed, especially two monolayer high pyramidal chains of atoms.In order to better understand the dynamics of nano structure formation, we present a kinetic Monte Carlo study on diffusion and adsorption of molybdenum atoms on a TiO₂(1 1 0) surface. A quasi one-dimensional lattice gas model has been used which describes the possible adsorption sites of a Mo atom in a single channel of the surface. The atomic positions of a 1.5 monolayer thick Mo film formed of pyramidal chains define the lattice sites of the model. Thereby the formation of three-dimensional clusters could be studied. Here we have studied the cluster formation as a function of parameters that can be controlled in a growth experiment by physical vapor deposition: deposition and annealing temperature, flux and total amount of deposited Mo. Good qualitative agreement with recent experiments is obtained.     

Deconvolution of very low primary energy SIMS depth profiles

In this paper, the deconvolution of SIMS profiles analysed at very low primary energy (0.5 keV/O₂⁺) is addressed. The depth resolution function (DRF) of the SIMS analysis in presence of roughness is established and a deconvolution procedure is implemented without or in presence of roughness on samples containing delta-doped layers of boron in silicon. It is shown that the deconvolution procedure can lead to a great improvement of the full width at half maximum (FWHM) of the measured peaks in the case where no roughness in detected in the crater bottom. In the case where it is present, the conditions required to use a deconvolution procedure are discussed, and the deconvolution is implemented using precise and restrictive assumptions.     

Au-analyte adducts resulting from single massive gold cluster impacts

Utilizing Au₄₀₀⁴⁺ primary ions produces large molecular ion yields, some in excess of unity, with minimal surface damage. A surprising observation is the occurrence of Au-analyte adducts as part of the ejecta desorbed by a single Au-cluster impact. We present data that demonstrate that Au and Au-adducts as secondary ions (e.g., AuCN⁻, AuGly⁻ and AuCsI⁻) are the result of the interaction between a single primary ion, Au₄₀₀⁴⁺ and the target atoms.     

Towards cooling of high-intensity ion beams

Buffer gas beam coolers may become excellent beam preparation devices for high-resolution mass separation. The small beam emittance provided makes efficient isobar resolution a realistic goal. In order to fulfill the needs of future facilities providing high-intensity beams of rare isotopes, it is desirable to increase the beam intensity limit of such devices from typically several tens of nanoamperes to microamperes. This requires the usage of high-voltage radiofrequencies in a low-pressure gas environment. A buffer gas beam cooler, dedicated to this purpose, is under development at the NSCL. The study of voltage breakdowns under such conditions and the design of an electrode system minimizing them is mandatory.      

Optimal vaccine scheduling in cancer immunotherapy

Cancer immunotherapy aims at stimulating the immune system to react against cancer stealth capabilities. It consists of repeatedly injecting small doses of a tumor-associated molecule one wants the immune system to recognize, until a consistent immune response directed against the tumor cells is observed.

We have applied the theory of optimal control to the problem of finding the optimal schedule of injections of an immunotherapeutic agent against cancer. The method employed works for a general ODE system and can be applied to find the optimal protocol in a variety of clinical problems where the kinetics of the drug or treatment and its influence on the normal physiologic functions have been described by a mathematical model.

We show that the choice of the cost function has dramatic effects on the kind of solution the optimization algorithm is able to find. This provides evidence that a careful ODE model and optimization schema must be designed by mathematicians and clinicians using their proper different perspectives.     

A systematic elimination procedure for Ito stochastic differential equations and the adiabatic approximation

We study a class of nonlinear Ito stochastic differential equations (with possibly state dependent diffusion coefficients), in which the variables can be divided into linearly damped (slaved) variables s and linearly undamped variablesu (order parameters). We devise a systematic and constructive procedure to eliminate the slaved variables. We take explicit time and chance dependence of the slaved variables into account, the latter via a family of diffusion processesZ _t ^(v) . These act as fluctuating coefficients of the Center Manifolds _t=s(u _t, t,Z _t ^(v) (v=2, 3, ...)) and appear explicitly in the elimination procedure. We show how in the Ito calculus fluctuating and deterministic coefficients of the Center Manifold are more completely separated than in the previously treated Stratonovich case [1]. The adiabatic approximation is defined as a partial summation of the elimination expansion and the stochastic generalization ofs=0 is derived. We show how thus ambiguity of stochastic calculi is removed. Closed form summations are given in two examples. We briefly indicate the potential use of perturbation theory techniques in the systematic elimination procedure.